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ABSTRACT 

We present a new method for constructing maps of the secondary temperature fluctu- 
ations imprinted on the cosmic microwave background (CMB) radiation by photons 
propagating through the evolving cosmic gravitational potential. Large cosmological 
N-body simulations are used to calculate the complete non-linear evolution of the 
peculiar gravitational potential. Tracing light rays back through the past lightcone 
of a chosen observer accurately captures the temperature perturbations generated by 
Hnear (the integrated Sachs- Wolfe or ISW effect) and non-linear (the Rees-Sciama or 
RS effect) evolution. These effects give rise to three kinds of non-linear features in the 
temperature maps, (a) In overdense regions, converging flows of matter induce cold 
spots of order 100 Mpc in extent which can dominate over the ISW effect at high 
redshift, and are surrounded by hot rings, (b) In underdense regions, the RS effect 
enhances ISW cold spots which can be surrounded by weak hot rings, (c) Transverse 
motions of large lumps of matter produce characteristic dipole features, consisting of 
adjacent hot and cold spots separated by a few tens of Megaparsecs. These non-Hnear 
features are not easily detectable; they modulate the ISW sky maps at about the 
10 per cent level. The RS effect causes the angular power spectrum to deviate from 
Hnear theory at / 50 and generates non-Gaussianity, skewing the one-point distribu- 
tion function to negative values. Cold spots of similar angular size, but much smaller 
amplitude than the CMB cold spot reported by Cruz et al. are produced. Joint analy- 
sis of our maps and the corresponding galaxy distribution may enable techniques to be 
developed to detect these non-linear, non-Gaussian features. Our maps are available 
at ,http: / / astro.dur.ac.uk/ ^cai/ISW , 

very large scales, k ^ 0.1 h Mpc~^ at 2; = 0, with the 
scale of this transition becoming even larger at higher red- 
shift. This discovery that the non- linear ef fect was m o re im - 
portant at higher reds hift was mad e_bv_^[Cai et al.1 (|2009l ) 
and later confirmed by ISmith et all ([2009). Recent studies 
using both perturbation theory and fitting formulae indi- 
cate that the combined ISW and RS effect makes a non- 
trivial contribution to the overall C MB non-G aussianit y as 
measured by the bispectrum (e.^. I Verde Sp ergel 20021: 
Giovi et al. I I2OO3I : IBoubekeur et alJl2009l : iMandlh Verdd 
2QQ9I ). The primordial-lensing-RS correlation contribution 
ma y yield an effective non -Gaussianity parameter, /nl, of 
10 (iMandlh Verdd 1 20091 ). Ah these results suggest that 
the RS effect is a very important supplement to the ISW 
effect even at very large scales. 

Observational evidence for a possible contribution from 
the RS effect at large scales has come from a combi ned study 
of th e SDSS LRG samples and the CMB (Gra nett et"al] 
I2QO9I V These authors report a 4-cr detection of the ISW 
and RS signal at the scale of 4 degrees, somewhat higher 
than expected from the standard flat ACDM universe. 
There have also been attempts to attribute an extreme 



1 INTRODUCTION 

The integrated Sachs- Wolfe (ISW) effect (ISachs Wolfd 
Il967l ) arises from the decay of the large-scale potential fluc- 
tuations as they are traversed by cosmic microwave back- 
ground (CMB) photons and induces secondary temperature 
perturbations in the CMB radiation. It occurs in both open 
cosmological models (Qm < 1) and models containing a cos- 
mological constant (Qa > 0) or dark energy. It provides an 
independent way of measuring the dynamical effect of dark 
energy, but it is challenging to detect: at large scales, the 
ISW signal suffers from cosmic variance, as there are too few 
independent modes, while at somewhat sm aller scales, it is 
entan gled with the Rees-Sciama (RS) effect (|Rees Sciamal 
ll968l V which arises from the non- linear evolution of gravi- 
tational potential perturbations. Since it is impossible to 
evade the cosmic variance at large scales, the information at 
small scales becomes valuable. To extract cosmological in- 
formation from such scales requires disentangling the ISW 
and RS effects. 

Studies using N-body simulations have found that the 
combined ISW plus RS power spectra of CMB tempera- 
ture fluctuations deviate from linear theory expectation at 
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cold spot in the CMB to the I SW and RS effects (e.g 
Martinez- Gonzalez Saii3 1990l: iMartinez-Gonzalez etaU 



1990l: iR udnick et al l 120071: llno ue Silk" 20061. l2007l : 



Tomita In oue 2008k iMasina Notari 2009a. b^. However, 



on theoretical grounds such an explanation seems unlikely 
because the estimated sizes of the non-linear structures re- 
sponsible for the cold spot are typically > 100 Mpc, which 



sumes Gaussian initial conditions (see alsolCruz et al.l 
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t al. 2005: Cruz et al.ll2006l: iMcEwen et al.1 
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I2OO7.: .McEwen et al.ll2008l. for discussion of non- 



Gaussianity.). However, it is still unclear whether or not the 
combined ISW and RS effects can generate cold spots of 
a few degrees with the right amplitudes and whether such 
large-scale non- Gaussianity can arise from large-scale struc- 
ture. A full understanding of the ISW and Rees-Sciama ef- 
fect could be crucial in explaining the oddities of these ob- 
servations. Meanwhile, the increase in sensitivity of forth- 
coming CMB experiments opens possibilities of exploit- 
ing CMB temperature fluctuations down to arcmin scales 
(Planclfl, ACtEI, SPtEI and APEX-S20), at which the ISW 
and RS effect may also entangle with other large-scale 
astrophysics of interest, i.e. lensing (e.g. I Verde Spergell 
l2002l : lNishizawa et ahllioosl : iMangilh Verde ll2009') and the 
Suny aev-Zerdovich (SZ) effect (ISunvaev ^sTz eldovich 1972) 
(e.g. ICooravll2002l : IPosalba et al.ll2003l : iBielbv et al.ll2009l ). 
To disentangle all these effects is the key to thoroughly ex- 
ploiting the information encoded in these upcoming CMB 
measurements. 

N-body simulations are the ideal tool for investigating 
the phenomena discussed above since they treat the non- 
linear regime accurately and permit the construction of a 
full sky map of the ISW and RS effects and the full un- 
derlying 3-dimensional light-cone. Maps of the ISW effect 
have be en constructed from both simulations a nd obser- 
vations (iBarreiro et al.ll2008l : ICranett et aLll2009l ). Most of 
them simply adopt the linear approximation, using only the 
density field to estimate the time derivative of the poten- 
tial. We will show in this paper that these linear maps are 
far from accurate. There are also maps constructed from 
ray-tracing through simulations (jTuluie Lagunal Il995l : 
IPuchades et al.ll2006l ). but these are limited by small simu- 
lation box sizes and are not adequate to explore very large- 
scale structures. A full sky map of the RS effect has been 
constructed using a constrained high-resolution hydrody- 
namical s imulation to model the RS effect in the very local 
universe (|Maturi et al.ll2007l V This map is useful for under- 
standing the RS effect from within the radial distance of 
110 Mpc, which is a very small volume. Maps from the ray- 
tracing of large cosmological volumes are still missing. 

In this paper, we develop a new method of constructing 
a full sky light-cone of the time derivative of the potential, 
using a large N-body simulation. Our method of com- 
puting 6 is fully non-linear and so should model the com- 
plete RS effect as well as the ISW component. Our Gpc box 
size simulation provides a sufficient number of independent 
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large scale modes to investigate the ISW effect fully. We 
ray-trace through the light-cone to produce maps of tem- 
perature fluctuations induced by the ISW and RS effects. 
Our maps cover a large range of scales and cosmic time 
with high accuracy and allow us to investigate the ISW and 
RS effects. The maps will also be a valuable source for un- 
derstanding CMB secondary non- Gaussianity arising from 
large-scale structure. They may also prove useful for disen- 
tangling lensing and SZ effects from the ISW and RS effects. 

The paper is organized as follows. In §2, we present the 
basic physics and the mathematical description of the ISW 
and RS effects. In §3, we describe our method of computing 
CMB temperature perturbations from our N-body simula- 
tion and ray-tracing to produce full sky maps from light- 
cone data. In §4, we identify and discuss three characteristic 
non-linear features of the temperature perturbations. Full 
sky maps are presented in §5. Finally, in §6, we discuss our 
results and draw conclusions. 



2 THE ISW AND REES-SCIAMA EFFECT 

In a ACDM universe, dominance of the cosmological con- 
stant, A, causes the expansion factor of the universe, a, to 
grow at a faster rate than the linear growth of density per- 
turbations, 5. Consequently, the cosmological constant has 
the direct dynamical effect of causing gravitational potential 
perturbations, $ oc —S/a, to decay. The ISW effect is caused 
by the change in energy of CMB photons as they traverse 
these linearly evolving potentials. A CMB photon passing 
through an overdense region, or cluster, will gain more en- 
ergy falling into the potential well than it later looses climb- 
ing out of the evolved shallower potential well. Therefore, 
overdense regions correspond to hot regions in a linear ISW 
map. The converse is true for a photon passing through an 
underdense region. Here, the potential fluctuation is posi- 
tive and the CMB photon looses more energy climbing the 
potential hill than it subsequently regains from its descent. 
Therefore, underdense regions appear cold in a linear ISW 
map. Non- linear growth of the density perturbations mod- 
ifies this picture, producing additional temperature pertur- 
bations - the RS effect. In overdense regions, the accelerated 
non-linear growth of structure acts to increase the depth of 
the potential wells resulting in a reduction in the CMB tem- 
perature, partially cancelling the ISW effect. In contrast, in 
underdense regions the RS effect enhances the ISW effect as 
saturation of the density contrast in voids further suppresses 
growth of the gravitational perturbation. We will discuss 
other situations in which the RS effect makes a significant 
contribution and analyze the morphology of the resulting 
features in the sky maps in § [S] 

The net induced ISW plus RS, hereafter ISWRS, tem- 
perature fluctuation along a direction h can be written as 
an integral of the time derivative of the gravitational poten- 
tial, j>, from the last scattering surface to the pres ent (i.e. 
ISachs Wolfelll967l : [Martinez-Gonzalez et al.lll990l ). 



AT(n) = -fc 



/ i(t. 



h) dt, 



(1) 



where t is cosmic time, tL the age of the universe at the last 
scattering surface, to the present age, i> the time derivative 
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of the gravitational potential, To the mean CMB temper- 
ature and c the speed of light. This is equivalent to the 
integral over radial comoving distance, r, 

AT(n) = ^fo f ^{r,n)a dr, (2) 



where is the comoving distance to the last scattering 
surface and a the expansion factor. 

To compute $ accurately from our simulation we make 
use of the Poisson equation expressed in comoving coordi- 
nates, \/^^{x,t) — 4:7TGp{t)a^6{x,t), which can be written 
in Fourier space as 



<i>{k,t) = 



Ho 
k 



5{k,t) 



(3) 



Here p{t) is the mean density of the universe, 6 the density 
contrast 6 = (p — p)/p, Hq and Qm the present values of 
the Hubble and matter density parameters and G is the 
gravitational constant. Taking the time derivative yields 



^sCk,t) 



S{k,t) 



(4) 



Combining this with the Fourier space form of the continuity 
equation, S{k, t) -\- ik ■ p{k, t) — 0, we have 



$(A:,t) 



k 



(5) 



where p{k,t) is the Fourier transform of the momentum 
density divided by the mean mass density, t) = [1 -\- 
S{x,t)]v{x,t). Equation (|5]) enables us to compute i> (in- 
cluding the contributions of both the linear ISW effect and 
the non-linear RS effects) to high accuracy using the density 
and momentum fields of our simulation. We will refer to the 
results obtained using equation (|5)) as the ISWRS. 

We wish to contrast these ISWRS predictions with 
the corresponding results from linear theory. In the linear 
regime, S{k,t) = D{t)6{k,z = 0), where D{t) is the linear 
growth factor. Substituting this into equation (|4| yields 

Hk,t) = I (^^y Qr.p{k,t)[l - mi (6) 

where /3{t) denotes the linear growth rate = 
d\n D(t)/d\na. This equation represents the conventional 
way of modelling the ISW effect and uses only the in- 
formation from the density field. It is equivalent to as- 
suming that the velocity field is related to the density 
field by the linear approximation p{k,t) — i5(k^t)k/k^ ~ 
i/3(t)S(k,t)(d/a)k/k'^ . We will refer to results obtained us- 
ing this linear approximation for the velocity field as LAV. 
Note that with this LAV approximation, $ simply scales 
with time according to the ISW linear growth factor, G{t) — 
dD{t)[l- /3{t)]/a^. 

In the LAV approximation the density field, directly 
determines the potential field, and its derivative, 
Hence the morphology of the $ field is determined directly 
by the form of the density field, 6, with overdense regions 
(S > 0) corresponding to positive regions of $ and under- 
dense regions corresponding to negative regions of In con- 
trast, for the exact calculation, represented by equation (|5l), 
this correspondence is broken and the dynamics of the den- 
sity field play an additional direct role in determining 
This is fully discussed in § 2] 



3 CONSTRUCTING FULL SKY MAPS 

In this section, we describe our method of constructing full 
sky maps of the ISWRS effect using our large simulation. 



3.1 The Gpc simulation 

To investigate the ISWRS effect, we need a simulation of suf- 
ficiently large volume to include the very large scale pertur- 
bation modes (hundreds of Megaparsecs) necessary to ver- 
ify convergence with linear theory. At the same time, we 
need sufficiently high resolution to investigate the effects of 
small-scale non-linearity on the structures resolved in the 
sky maps. Moreover, for high accuracy ray- tracing, we need 
the redshift spacing of the simulation outputs to be small 
enough that interpolation between neighbouring redshifts 
does not introduce significant systematic errors. 

The N-body simulation we employ follows 2200^- 
particles in a 1 h~^Gpc periodic box in a ACDM cosmologi- 
cal model with Qa = 0.74, Qm = 0.26, Qb = 0.044, ag = 0.8 
and Ho = 71.5 km s~^ Mpc~^, chosen for co nsistency with 
recen t CMB and large scale structure data (| Sanchez et al.l 
l2009l ). The simulation, run on the COSMA supercomputer 
at Durham, was designed in order to make mock galaxy 
catalogues for forthcoming surveys (e.g. Pan-STARRSl and 
EUCLID) and, as such, resolves the dark matter halos of lu- 
minous galaxies (Baugh et al. in prep). The simulation has 
a softening length (Plummer equivalent) of 0.023 h~^Mpc 
which provides more than adequate resolution for our pur- 
poses. The initial conditions were set up at redshift z — 49 
and the simulation was run using GADGET (|Springel et al.l 
[2005) with output at 50 snapshots, 48 of which lay between 
z — and z — 10, where we mainly focus our analysis. 
The redshift intervals between neighbouring simulation out- 
puts correspond to about 100 Mpc in radial comoving 
distance. We have verified that this is adequate to ensure 
that the errors induced by interpolating i> between snap- 
shots are less than 2% at Megaparsec scales and even less 
at larger scales. 



3.2 Map construction 

Our method of constructing the ISWRS sky maps has two 
stages. First, at each output redshift, we construct an esti- 
mate of i> on a cubic grid. In the second stage, we propagate 
light rays from the observer and, as we move along each ray, 
we interpolate $ from the grids and accumulate the AT 
defined by the integral of equation (|2]) along t he past light- 
cone of the observer. Finally, we use HEALPix (|G6rski et al.l 
2005) to visualize the AT map in spherical coordinates. We 
now describe these steps in more detail. 



3.2.1 Cartesian Grids 

We construct the density field, ^(x), by assigning the 
mass of each dark matter particle to a 3D mesh of cu- 
bic grid cells using the cloud-in-cell assignment scheme 
(|Hocknev E astwood 1981). We apply the same scheme to 
accumulate the momentum density field, p(x), by assigning 
the vector momentum of each particle to a Cartesian grid. 
Then we perform four independent Fast Fourier Transforms 
to compute the Fourier space versions of the density field. 
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S{k), and the three components of the momentum field, p{k). 
These are then combined using equation ([5|) to yield the 6 
field in Fourier space. Finally, we perform an inverse Fourier 
transform to obtain $(x) in real space on a cubic grid. We 
repeat the final two operations using equation (|6|) to ob- 
tain the alternative LAV $ field. This whole process is then 
repeated for each of the 50 simulation outputs. 

For our 1 Gpc simulation, we have used a grid of 
1000"^ cells each of 1 Mpc on a side. The resolution of the 
N-body simulation is significantly greater and would warrant 
using a finer grid. However, this would be very demanding 
of both machine memory and hard disk space and is not 
necessary to resolve structures over 10 Mpc accurately, 
which is the smallest scale of interest in this paper. 



3.2.2 Ray tracing through the light-cone 

The next step is to choose a location for the observer and 
propagate rays through the simulation. We choose to place 
the observer at the corner of the simulation bo x at Cartesian 
coordinate (0,0,0). We then used HEALPix (iGorski et al.1 
I2OO5I I to generate the directions of 3 145 728 evenly dis- 
tributed rays, corresponding to an angular pixel scale of 
(6.87^)^ which is sufficient to resolve upto a spherical har- 
monic scale of / ^ 1000. For each ray, we accumulate the 
integral given by equation ([2|) by taking fixed discrete steps 
in comoving radial distance. At the location of each step, 
we find the two output snapshots that bracket the lookback 
time at this distance. Using the Cartesian grids at each of 
these outputs we use the c loud-in-cell assignment scheme 
(iHocknev Eastwoodlll98lh to obtain the values of $1 and 
<i>2 at the chosen position on the ray at these two lookback 
times. Note that if the position along the ray lies outside the 
simulation box, we use the periodic boundary conditions to 
map the location back into the box. Finally, to estimate $ 
at the lookback time corresponding to the position on the 
ray, we linearly interpolate ^{t)/G{t) with comoving radial 
distance, r, using 

(r2 - ^i)§ = (^2 - r)^ + (r - n)^. (7) 

Here C, Gi and G2 are the linear growth factors for i> at 
the lookback time corresponding to the position on the light 
ray and the two neighbouring outputs that bracket it. The 
values of r, ri and r2 are the corresponding comoving radial 
distances. This interpolation scheme guarantees that we re- 
cover the linear theory result exactly when $ is evolving 
according to linear theory. 

The maps that we present in Section [5] ar e constructed 
and analysed using the HEALPix package (jCorski et al.l 
I2OO5I I. We show maps corresponding to the contribution of 
the integral in equation (|2]) over different finite intervals of 
comoving radial distance. If we were to integrate over a co- 
moving distance larger than the 1 h~^Gpc size of our simula- 
tion, the assumed periodic boundary conditions would create 
artifacts in the maps. For instance, in directions correspond- 
ing to the principal axes of the simulation volume the light 
ray would pass through the same location every 1 /i~^Gpc. 
The contributions to AT(n) from these replicas would add 
coherently and lead to larger fluctuations along these spe- 
cial directions than on average. In similar applications some 



authors choose to rotate, flip and shift the replicated simula- 
tion cube. While this avoids the artificial coherence between 
one replica and the next, it generates other problems includ- 
ing discontinuities of i> at cube boundaries. We have chosen 
to evade this difficulty by mainly focusing on generating and 
analyzing maps with a maximum radial depth equal to the 
simulation cube size of 1 h-^Gpc. It may still remain the 
case that the same structure is seen in more than one di- 
rection, but this does not affect the mean power spectrum 
of the map and only produces non-Gaussian features on an- 
gular scales approaching the angular size subtended by the 
simulation cube. Hence, the power spectra and small scale 
features of the maps are not affected by the periodic nature 
of the simulationjj 

It is also important to establish that the resolution of 
the maps is not compromised by the spatial and temporal 
resolution of the simulation or our choices of size of grid cell 
and integration step. The pixel size of our HEALpix maps 
is 6.87^ and matches the linear size of our Cartesian grid 
cells at an angular diameter distance of 500 Mpc. Thus, 
we would expect that beyond 500 Mpc our sky maps 
are not lacking any small scale features due to the limited 
resolution of our cubical grids. We have tested this using a 
smaller simulation of the same resolution but grided using 
0.5 rather than 1 Mpc cells. This test indicates that at 
/ = 100, the larger, default, choice of cell size suppresses the 
power in the maps made from the innermost 250 h~^Mpc 
by just 10%. As we ray-trace further in the radial direc- 
tion, the accuracy improves, and is within 10% for / ^ 1000 
for maps extending to 1000 h~^Mpc. Additional tests have 
shown that it is this cell size which is the limiting factor in 
the resolution of our maps. For instance, we have checked 
that the maps are essentially unchanged if we make the in- 
tegration step size smaller than our default choice or if we 
have more closely spaced simulation outputs. The full reso- 
lution of the N-body simulation could be exploited by using 
a finer mesh but this is unnecessary for our purposes. 



4 THE EFFECTS OF NON-LINEARITY 

In this section, we compare our full non-linear estimates of 
the ISWRS effect with those of the LAV approximation. We 
both quantify and characterise the features that are gener- 
ated by non-linear gravitational evolution and elucidate the 
physical processes by which they are generated. 

Fig. [1] shows the temperature maps at three different 
redshifts that result from computing the contribution to the 
ISWRS integral equation ^ from a 100/i~^ Mpc thick shce 
of the Gpc simulation. The top row shows the full ISWRS 
calculation in which 4> is computed using equation (|5)) , while 
the middle row shows the result of using the LAV approxi- 
mation for $ given by equation (|6| . The LAV differs slightly 



^ In fact, we find that if we ignore this problem and produce 
power spectra directly from maps projected over a depth of 
7000 h~^Mpc or more, they do not differ significantly from those 
built up from the power spectra of successive 1000 h~^Mpc slices. 
Hence, at least for angular power spectra, our simulation size is 
sufficiently large and the evolution of <!> sufficiently rapid that the 
superposition of periodic replicas is not a concern. 
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Figure 1. Maps of AT generated from a slab of thickness Ar = 100/i~^ Mpc taken from our Gpc simulation. From left to right, they are 
maps at z = 0.0, 2.1, 6.2 respectively. From top to bottom, they are maps constructed using equation (O, which includes the ISW and 
Rees-Sciama effects (ISWRS), maps constructed using the linear approximation for the velocity field, equation (|6]), (LAV), and residual 
maps of the top panels minus the middle panels, leaving essentially the Rees-Sciama (RS) contribution to the temperature fiuctuations. 
Note the individual temperature scales for each panel. At z = 0, we also show the momentum field, averaged over the same slice of the 
simulation, by the overplotted arrows. The three square boxes indicate the three regions we show in detail in Fig. [3] (solid-black box). 
Fig. [5] (dotted-black box) and Fig. [6] (dashed-yellow box). 



from the standard ISW contribution to the temperature fluc- 
tuations because to compute the ISW one should really use 
the linear theory prediction for the density field in equa- 
tion (|6|), whereas for the LAV we use the actual non-linear 
density field, but assume the velocity field is related to the 
density as in linear theory. However, the difference between 
the LAV and the true ISW contribution is extremely small, 



as we shall see when we compare their power spectra in Sec- 
tion \5\ Consequently, the bottom row of panels, which are 
the residual produced by subtracting the LAV maps from 
the corresponding ISWRS maps, are essentially the RS con- 
tribution to the temperature fluctuations. 

Let us consider first the LAV maps. It is clear that the 
hot (red/yellow) and cold (blue/green) regions have a large 
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Figure 2. Temperature perturbations along three light rays, 
shown in black, blue and red from 2; = 0to2; = 5.8, correspond- 
ing to the comoving distance from Tc = to Tc = 6000 Mpc. 
The top three panels show the temperature perturbations per 
unit comoving distance along each radial direction. Solid lines 
show the results of the full ISW and Rees-Sciama effects (ISWRS) 
and dashed lines the results of the linear approximation for the 
velocity field (LAV). The fourth panel from the top shows the ac- 
cumulated temperature perturbations along the radial direction 
for ISWRS (solid lines) and LAV (dash lines). The bottom panel 
shows differences between these accumulated ISWRS and LAV 
temperature perturbations. 



coherence scale, approaching the size of the simulation box. 
The coherence scale is much larger than that of the un- 
derlying density field simply because of the factor in 
equation (|4|, which arises from the Poisson equation relat- 
ing $ and 6. As explained at the beginning of Section O 
hot and cold spots in these maps are induced by the de- 
cay of the gravitational potential in overdense and under- 
dense regions respectively caused by the dynamical effect of 
the cosmological constant, Qa- In the higher redshift slices, 
the amplitude of the LAV/IS W fluctuations are greatly re- 
duced. This is easily understood: as redshift increases Qa 
becomes smaller and linear theory predicts density pertur- 
bations grow as ^ oc a and so the corresponding potential 
perturbations, $ oc —6 /a, are constant. For 6 = there 



are no ISW fluctuations as CMB photons traversing such 
a potential well will loose just as much energy leaving the 
perturbation as they gain on entering it. 

The RS contribution to the temperature fluctuations, 
shown in the bottom row of Fig. [U have a much shorter 
coherence scale than the ISW fluctuations. Their evolution 
with redshift is much more gradual than for the ISW fluc- 
tuations. In the redshift z = slice they are a minor con- 
tribution to the overall ISWRS map (top row), but they 
become a significant contribution at redshift z = 2.1, gen- 
erating small scale structure within the smooth large-scale 
ISW fluctuations. At z — 6.2, the RS fluctuations are al- 
most completely dominant and while they are on a smaller 
scale than the ISW fluctuations, they still produce filamen- 
tary structures which can be seen to be one hundred to a 
few hundreds of Megaparsecs across. These panels give a 
visual confirmat ion of the two-point statistics presented in 
ICai et al.1 (l2009l ). which demonstrated that both the impor- 
tance and physical scale of the RS effect become greater with 
increasing redshift. 

Fig. [2] compares the ISWRS and LAV temperature per- 
turbations along three randomly chosen light rays (avoiding 
the principal axes of the simulation box) . In each of the top 
three panels the solid line shows the overall ISWRS pertur- 
bation and the smoother dashed line the LAV contribution 
from z = 0to^ = 5.8, corresponding to the comoving dis- 
tance from Tc = to Tc = 6000 Mpc. As in the maps 
of Fig. [1] we see that the LAV contribution varies more 
smoothly and rapidly damps in amplitude as one progresses 
along the rays to higher redshift. The RS contribution pro- 
duces high frequency fluctuations around the LAV predic- 
tions with the amplitude of these high frequency fluctuations 
decreasing with redshift at a much more modest rate. The 
lower two panels of Fig. [2] show that while a large portion of 
the high frequency RS signal cancels out when accumulated 
along the line- of- sight, persistent residuals of the order of 
1 /iK are also accumulated with a large contribution coming 
from the redshift range around a radial distance of about 
2000 h-^ Mpc. 

In order to identify the non-linear physical processes 
that give rise to the persistent features in the ISWRS maps, 
we have studied the momentum field of our simulation. The 
overlaid white arrows in the left hand panels of Fig. [1] indi- 
cate the projected momentum field in this slice of the sim- 
ulation. At higher redshifts we find that the orientation of 
the momentum vectors is similar to that at 2; = 0, but their 
amplitude increases, roughly in accord with the expected 
linear growth rate factor, {3. In the LAV panel (middle row) 
we see that on large scales the momentum field is well cor- 
related with the temperature perturbations and hence with 
the $ field. Dark matter particles are moving towards hot 
lumps (overdense regions, potential wells) and flowing away 
from cold lumps (underdense regions, potential hills). This 
is the scenario expected in linear theory. However, looking at 
the smaller scale features evident in the ISWRS (top panel) 
maps we flnd that this correlation between temperature and 
momentum does not always hold. In fact, at high redshift, 
violations of the correspondence expected in linear theory 
are very strong. Hence, guided by the momentum fleld we 
have identifled three interesting non-linear phenomena in the 
ISWRS maps, which are related to dipoles, convergent flows 
and divergent flows. Three regions exhibiting these phenom- 
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Figure 3. Maps of AT in a slice of thickness Ar^SO/i"-*^ Mpc and an area of 50 x 50 \h~^ Mpc]^. The x and y labels give the exact 
coordinates of this set of maps in the full-box maps shown in Fig. [T] From left to right, they are maps at z = 0.0, 2.1, 4.2 respectively. 
From the top to the bottom, they are maps constructed using equation (|5]), which includes the ISW and Rees-Sciama effect (ISWRS), 
maps with linear approximation for the velocity field, equation (|6]), (LAV), and residual maps of the top panels minus the middle panels, 
which is essentially the Rees-Sciama (RS) contribution. The overplotted arrows show the projected momentum field of the slice. 



ena are highlighted by the square boxes in Fig. [T] (one of 
which is split across the periodic boundary of the simula- 
tion). We zoom in and study each in detail in the following 
subsections. 

4.1 Dipoles 

In RS maps, at scales of tens of Megaparsecs, large lumps of 
dark matter moving perpendicular to the line-of-sight give 



rise to dipole features, i.e. a cold spot on the leading part of 
the lump and a hot spot on the trailing part. Fig. [3] shows 
zoomed-in maps of the ISWRS, LAV and RS temperature 
perturbations in the region around one such dipole at red- 
shifts ^ = 0, 2.1 and 4.2 for a 50 /i"^Mpc. The overplotted 
momentum vectors clearly show the bulk motion of a large 
lump of dark matter. The length of the arrows indicate that 
this velocity is being damped and so decreases with decreas- 
ing redshift. The LAV maps are very smooth and show no 
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Figure 4. Segments of the three hght rays shown in Fig.O with the distances shown on the x-axis. The top panels show the temperature 
perturbations per unit comoving distance in these chosen directions. Sohd hues show the results of the full ISW and Rees-Sciama effects 
(ISWRS) and dashed lines those from the linear approximation for the velocity field (LAV). The middle panels show the accumulated 
temperature perturbations along the radial direction for ISWRS (solid lines) and LAV (dashed lines). The bottom panels shows the 
accumulated difference between the ISWRS and LAV perturbation. The left panel is an example of a dipole, the middle panel of a 
convergent flow and the right panel of a divergent flow around an empty void. 



sign of the dipole feature at all, while the strikingly large 
dipole with amplitude of 1 /iK is clearly visible in the 
ISWRS maps (top row), particularly at the higher redshifts. 

The physical origin of these dipole features, which are 
just a special case of the RS effect, can be understood in 
terms of the evolution of the gravitational potential, on 
either side of the moving mass. At a fixed position ahead 
of the moving mass, $ decreases as the mass and its gravi- 
tational potential well approaches. This will create a CMB 
cold spot as CMB photons passing through this point will 
gain less energy falling into the potential well than they 
will subsequently loose climbing out of the then deeper well. 
Conversely, behind the mass, $ is increasing (becoming less 
deep) and so CMB photons gain more energy than they loose 
and this creates a hot spot. 

As the maps in Fig. [3] are just the contributions to 
AT from a thin, 50 /i~^Mpc, slice of our simulation, one 
might worry that this artificial truncation gives rise to ar- 
tificial edge effects. We have checked that this is not the 
case by making corresponding maps projected through a 
depth of 1 /i~^Gpc, the box size of the periodic simula- 
tion. We find that the dipole features are still clearly vis- 
ible, but, of course, they are somewhat perturbed by su- 
perposition of other perturbations along the longer line-of- 
sight. This finding is consistent with the accumulated dif- 
ference between the ISWRS and LAV AT shown earlier in 
Fig- [2] (bottom panel), which remain roughly constant for 
Tc > 3000 /i~^Mpc in spite of some small-scale variations. 

The amplitude of the temperature fluctuations 



generated by the dipoles is not large, but their char- 
acteristic dipole sig nature ni i ght e nable them to be 
detected (.Rubino-Martm et al.1 l2004l : iMaturi et all l2006l . 
l2007l ). Dipoles generated by moving galaxy clusters, 
which can equally be thought of a s ^mov ing lenses', 
were predicted bv iBirkinshaw Gulll (Il983h , and later 
discussed by ICurvits Mitrofanovl (|l986h . Their de- 
tect ability in the CMB and possible uses, for example to 
measure transverse motions of dark matt er , have been 
further exp lored bv [ Tuluie Lagun 1 ("1995"); Tuluie et al. 
(199g); .Aghanim et al.l (Il998l); [M olnar & Birki nshaw 
(2 00(1: lAso et al.1 ([2002lV . 
iRubifio-Martm et al.1 (^Ooj ) 



Molnar Birkinshawl 
ICoorav Setd 
iMaturi et al.l (120071 ). These studies considered only dipoles 
on the scale of galaxy clusters while our simulations 
reveal dipoles of /jK amplitude on scales ranging from 
10 h~^Mpc to a few tens of /i~^Mpc, the larger of which 
are seeded by bulk motions on scales far larger than galaxy 
clusters. 

Moving dark matter lumps will perturb the energies of 
CMB photons even if they move along, rather than trans- 
verse to the line- of- sight. Examples of the perturbations such 
bulk motions cause are evident at several points along the 
rays shown in Fig. [2] They appear as a sharp peak followed 
very closely by a deep dip. The left-hand panels of Fig. |4] 
show a close-up view of one of these features. Following the 
ray on the top panel from the right to the left (i.e. moving 
in the direction of a CMB photon), the line dips and then 
peaks, indicating that a lump of dark matter is moving in 
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Figure 5. Like Fig. O but for slice of thickness Ar=100h~^ Mpc, with an area 
From left to right the panels show maps at z = 0.0, 2.1, 6.2 respectively. Note 
simulation box (see Fig. To show it properly, we shifted the simulation box 
boundary conditions, so that the y-axis ranges shown in the maps correspond 
90h~^ Mpc in the original simulation cube. 
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of 200 X 200 [h~^ Mpc]^ centred on a convergent flow, 
that the region shown here lies at the boundary of the 
along the y-axis by —100 Mpc using the periodic 
to the combination of 890 to 1000/i~^ Mpc and to 



the opposite direction to the CMB photon. The local po- 
tential on the lump's leading part is getting deeper (cooling 
down photons) while on the trailing edge the potential is 
becoming shallower (heating up photons). The accumulated 
AT (middle panel) is boosted and then suppressed as the ray 
passes through and the net effect is extremely small. This is 
to be contrasted to the case discussed previously when the 
lump moves transverse to the line-of-sight for which there is 
no cancellation. 



4.2 Convergent Flows 

In the RS and the higher redshift ISWRS maps of Fig. □ 
we can see several examples of cold regions surrounded 
by hot rings or filaments which are centred on converg- 
ing velocity flows. These features are larger in scale than 
the dipoles, ranging from tens to hundreds of Megaparsecs. 
Fig. [5] shows zoomed-in maps of the ISWRS, LAV and RS 
temperature perturbations and overlaid velocity vectors at 
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Figure 6. As Fig. [3 but for a slice of thickness Ar=1000h-^ Mpc, with an area of 300 x 300 [h' 
From left to right the panels show maps at z = 0.0, 2.1, 6.2 respectively. 
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Mpc]^ centred on a divergent flow. 



redshifts z = 0, 2.1 and 6.2 for a 200 /i~^Mpc box centred 
on one such feature selected from Fig. [1] In the LAV maps 
(middle row), we see only a smooth hot lump, centred on 
the convergent flow, whose amplitude increases towards low 
redshift as Qa becomes increasingly important. In contrast, 
the RS contribution to the temperature perturbations (bot- 
tom row) is negative at the centre of the convergent flow and 
is surrounded by a positive shell. Its amplitude evolves with 
redshift much more weakly than the amplitude of the LAV 
perturbation. It is strongest in the central panel, decreas- 
ing at the lowest redshift due to the damping of non-linear 



growth by the late time acceleration of the universe. At all 
but the lowest redshift, this RS perturbation is a significant 
contribution to the overall ISWRS maps (top row). At high 
redshift, the RS completely dominates over the LAV contri- 
bution turning the hot spot predicted by linear theory into 
a cold spot, almost 200 /i~^Mpc in extent, with amplitude 
of order a /xK, surrounded by a hot shell. Even in the red- 
shift z = 2.1 slice, the RS contribution drastically changes 
the morphology of the LAV hotspot, producing a cold re- 
gion in the very centre of the fiow that is surrounded a hot 
filamentary shell. 
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The explanation of this counter-intuitive conclusion 
that overdense regions can become cold spots surrounded by 
hot-rings, rather than the hotspots predicted by the ISW, 
involves principally the same physics as is responsible for 
the dipoles. In the ACDM model, overdense regions grow 
as the result of the inflow of lumpy material, often along 
filaments. Each of these inflowing lumps will give rise to a 
dipole feature. On the leading edge of the lump, the poten- 
tial is decreasing and CMB photons loose energy, while on 
the trailing edge the potential is increasing and CMB pho- 
tons are heated. The only difference is that dipoles seen on 
the sky are indicators of lumps of dark matter with large 
transverse momentum, while these larger scale cold regions 
surrounded by hot rings are produced by larger scale coher- 
ent convergent flows. One can imagine splitting the conver- 
gence flow into many small lumps of dark matter moving to- 
wards the same centre: the cold region in the centre is just 
the result of stacking many leading parts of those lumps, 
while the hot ring consists of their trailing parts. Our find- 
ings regarding the morphology of the RS effect in converg- 
ing flows are in general agreement with analytic models of 
form ing clusters that have been discu s sed by other authors 
e.g. [Martinez- Gonzalez Sand 1990l: iLasenbv et aPllQQOl : 



iDabrowski et al.lll999l : llnoue Silkll2006l . l2007l ). 

The maps shown in Fig. [5] are the contribution to the 
AT perturbation of just a 100 h~^Mpc slice of our simu- 
lation. Again, one should consider whether this truncation 
gives rise to artificial edge effects that would qualitatively 
affect the appearance of the cold spots and hot rings. In 
particular, if the slice artificially removes a foreground or 
background section of the hot shell this will enhance the 
visibility of the cold central region. We have directly tested 
whether this is a strong effect by making untruncated maps 
whose depth is the full 1000 h~^Mpc size of the periodic 
simulation. We find that while there are some cases where 
the heating and cooling cancel each other out, normally the 
visibility and contrast of the cold spot and hot ring features 
is not strongly affected. This is perhaps to be expected as 
the hot shell is diluted because it is spread over a much 
larger area than is subtended by the cold spot. Moreover, 
the dominance of the cold spots can be seen directly in the 
plots of the AT contributions along particular lines of sight. 
An example of a ray passing through a cold spot is shown 
by the middle column of Fig. [U In the top panel the broad 
peak predicted by the LAV (dashed line) corresponds to an 
overdense region in the simulation. There is a convergent 
flow around this overdense region that gives rise to both 
the sharp central dip and surrounding upward fluctuations 
seen in the ISWRS result (solid line). The lower panels show 
that the difference between ISWRS and LAV perturbations 
is strongly dominated by the central cold spot. 

4.3 Divergent Flows 

Divergent flows surrounding voids or underdense regions also 
produce characteristic features in the ISWRS and RS maps. 
However their effect is not simply the reverse of that of 
the convergent flows. Instead, the non-linear behaviour in 
these void regions always acts to enhance the LAV perturba- 
tion producing stronger cold spots. Fig. [6] shows zoomed-in 
ISWRS, LAV and RS temperature maps of a 300 h-^Mpc 
region centred on such a cold spot. This region is underdense 



and the LAV maps (middle row) show the expected linear 
ISW behaviour of a cold spot which grows in amplitude with 
decreasing redshift as Qa increases and the potential hill as- 
sociated with the underdensity decays. The RS contribution 
(bottom row) also produces a cold spot at the centre of the 
region, but surrounded by a hot filamentary shell. Thus, per- 
haps counter-intuitively, the pattern of the RS contribution 
for these divergent flows is the same, and has the same sign, 
as for the convergent flows discussed above. 

The reason why the effect of non-linearity in an under- 
dense region is not simply the reverse of its effect in an over- 
dense region is because the two situations are not symmet- 
rical. In an overdense region the overdensity is unbounded 
while in a underdense region the density cannot drop below 
zero, so 5 > —1. It is this saturation of the density contrast 
in voids that prevents the perturbation from growing in the 
non-linear regime even as fast as linear theory would predict 
and so enhances the rate of decay of 

Again, the amplitude of the RS contribution evolves 
with redshift much more weakly than the LAV contribu- 
tion. At 2; = 0, the LAV contribution to the ISWRS is over- 
whelmingly dominant, but at higher redshift, the RS effect 
contributes by reinforcing the linear effect in the centre of 
the cold spot and suppressing it in the outer regions. 

4.4 Overview of non-linearity 

A useful summary of the three physical effects we have dis- 
cussed above is given by the schematic diagrams in Fig. [71 
They depict the gravitational potential, $ oc —S/a, at a 
baseline reference epoch, say when the CMB photon enters 
the region, and both the LAV and fully non-linear predic- 
tions for the evolved potential at a later time, say when 
the CMB photon exits the region. These schematics yield 
straightforward interpretations of each of the phenomena 
we have identified in the ISWRS and RS maps. 

In void regions, the ISW effect is easily understood as 
the result of the linear decay of the gravitational potential. 
CMB photons are cooled as they loose more energy climbing 
the initial potential hill of the void region than they subse- 
quently regain on departing the now shallower potential hill. 
This cooling is depicted by the heavy arrows in the left-hand 
panel of Fig. [71 One would expect non-linearity to acceler- 
ate the growth of 6 and hence increase the rate of growth 
of the potential hill. This is what we observe to occur at 
the edge of the voids and leads to a component of heating 
of the CMB photons in this outer shell, as depicted by the 
light upward arrows. However, once the centre of a void re- 
gion becomes almost empty, i.e. ^ ~ — 1, the local density 
contrast stops growing, as it cannot become emptier. In this 
case, the expansion of the universe will reduce the height of 
the potential hill, just as it does in the linear regime in voids. 
Thus, in the centre of voids the RS effect has the same sign 
as the linear ISW effect and the two reinforce each other to 
produce cold spots at the centre of voids. 

In overdense regions such as galaxy clusters, the lin- 
ear decay of the gravitational potential well results in the 
heating of CMB photons as depicted by the heavy arrows 
in the central panel of Fig. [71 The effect of non-linear in- 
fall and growth of the density perturbation is to shrink the 
scale of the potential well while deepening its central value. 
In the centre of the cluster this non-linear growth acts in 



12 Cai et al. 




Figure 7. Schematic diagrams of the evolution of the gravitational potentials that produce the three characteristic features of the 
non-linear RS effect discussed in Section (4] The solid lines depict for a void (left), cluster (middle) and moving cluster (right) the 
gravitational potential at a reference epoch, say when a CMB photon enters the potential. The dashed and dotted lines show predictions 
for the evolved potential at a later time, say when the CMB photon exits the region. The dashed lines are the linear prediction of the 
LAV approximation and the dotted line represents the fully non-linear result. The heavy up and down arrows indicate the heating or 
cooling of the CMB that is predicted by the LAV approximation (essentially the linear ISW effect). The lighter arrows indicate the RS 
contribution to the temperature perturbations, i.e. the difference between the fully non-linear and LAV contributions. 



opposition to the linear decay of the potential and cools the 
CMB photons, depicted by the downward light arrow. We 
have seen that at high redshift this effect can be stronger 
than the linear ISW effect, resulting in cold spots centred 
on overdense regions. At larger scales, the inflowing region 
shrinks the scale of the potential well and results in a further 
reduction in the depth of the potential well, reinforcing the 
linear effect. Hence, in an outer shell around the cluster, the 
non-linear effect heats CMB photons. This latter effect can 
also be considered as the superposition of a set of dipoles 
arising from lumps in a surrounding inflowing region. 

The right-hand panel of Fig. [71 depicts the situation that 
gives rise to a characteristic dipole perturbation of the CMB. 
A moving cluster, or other overdense region, gives rise to a 
moving potential well. Linear theory does not model the 
movement of the potential well and would instead predict a 
simple static decaying potential well, resulting in a hot spot 
as indicated by the heavy upward arrows. In contrast, the 
movement of the potential well leads to a rapid deepening 
of the potential well on the leading edge of the cluster and a 
corresponding rapid increase of the potential on the trailing 
edge. This gives rise to cooling of the CMB photons on the 
leading edge and heating on the trailing edge as indicated 
by the light arrows. 



5 SKY MAPS 

We have used the methods described in 33] to construct 
full sky maps of the ISWRS effect using both the density 
and velocity fields, and using the LAV approximation. We 
have also generated maps of the RS effect, by subtracting 
the LAV maps from the corresponding ISWRS maps. In 
all of the sky maps we have removed both the monopole 
and dipole using the REMO VE_DIPOLE subroutine of 
HEALPix (jGorski et al.ll2005l ). 



The top panel of Fig. [8] shows the ISWRS map made 
by integrating along the line-of-sight of the observer over 
the range < z < 0.17, corresponding to distances from the 
observer, Tc, in the range < rc <500/i~^ Mpc. The whole 
sky is dominated by a few large hot and cold features with 
amplitudes of a few /iK to 10 /iK, as expected from linear 
theory. The bottom panel shows the ISWRS map integrated 
from 0.17 < z < 0.57, corresponding to 500 < rc < 1500 

Mpc. The typical angular size of the features in this 
shell is smaller than in the top panel. As we will discuss later 
in this subsection, the smaller angular size of the features 
in the bottom map relative to the top map is due largely 
to the fact that there is a cut-off in the power spectrum of 
the simulation used to build the sky maps. Unlike in the 
real universe, there is no power contributing to our maps on 
scales larger than the fundamental modes of the simulation. 
This defect can easily be remedied: as the missing modes are 
essentially in the linear regime, it would be straightforward 
simply to add additional longer wavelength modes by hand 
before doing the line-of-sight projection to make the maps. 
This extra power, however, would be the same in maps of 
the ISWRS effect and the LAV approximation, and so would 
vanish identically from the maps of the RS effect. For this 
reason we have chosen not to add these linear modes to the 
ISWRS and LAV maps. 

The corresponding maps, but with the LAV approxima- 
tion, are shown in Fig. [9l Comparing Fig. [8] with Fig. [9l we 
see that the large-scale distribution is essentially identical, 
but differences are apparent on smaller scales. The difference 
map made by subtracting the LAV map from the ISWRS 
map to reveal the RS effect, is shown in Fig. 1101 In the top 
panel, covering the range < 2; < 0.17, we see some striking 
large-scale structures with amplitudes that are about l/xK, 
which is about 10% of the amplitude of the features seen in 
the ISWRS and LAV maps themselves. A few large dipoles, 
ranging from a few degrees to over 10 degrees are clearly 
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Figure 8. Full sky maps of the predicted CMB AT due to the ISW and Rees-Sciama effects (ISWRS) made from our simulation. The 
upper plot is obtained by ray-tracing through the simulation over the redshift interval < z < 0.17, which corresponds to a range of 
comoving distance from the observer of - 500 Mpc. In the lower plot the projection is for the range 0.17 < z < 0.57 (corresponding 
to the range 500 - 1500 Mpc in comoving distance). These maps, and all subsequent sky maps, use the Mollweide projection to 
represent the sky on a plane, with each pixel having an area of (6.87^)^. The grid spacing is 30° in both longitude and latitude. For 
reference, the box-size of our simulation is 1000 Mdc. 
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Figure 9. As Fig. [S] but showing the full sky maps of the predicted CMB AT due to the ISW effect constructed using the linear 
approximation for the velocity field (LAV). 
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Figure 10. As Fig. [S] but for the Rees-Sciama (RS) effect, obtained by subtracting the maps in Fig. [9] from the corresponding maps in 
Fig. El 
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Figure 11. Maps of 40 x 40 degree patches of the sky of, from left to right: (i) the ISW and Rees-Sciama effects (ISWRS); (ii) the 
hnear approximation for the velocity field (LAV); (iii) the Rees-Sciama (RS) effect. The projection is over the range 0.57 < z < 1.07, 
corresponding to 1500 - 2500 Mpc comoving distance from the observer. The grid spacing is 10° in both longitude and latitude. 
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Figure 12. As Fig.[TT] but projecting over the range 1.07 < z < 1.78, corresponding to a comoving distance range of 2500 - 3500 Mpc. 
Note that since the simulation volume is only 1000 Mpc on a side, the same features necessarily appear, but at smaller angular 
scales than in the previous figure. 



visible, indicating the large-scale bulk flow of matter in the 
local universe of our observer. The bottom panel shows the 
RS effect for a projection covering the range 0.17 < z < 0.57. 
Many more dipoles are visible, with typical angular scales of 
a few degrees. In addition, the map has large-scale filamen- 
tary structure which is coherent over lengths of up to tens 
of degrees. The typical amplitude of these features is a few 
/iK, again about 10% of the amplitude of the features in the 
corresponding ISWRS and LAV maps. 

In Fig. [11] and Fig. [121 we show a patch of sky of size 
40 X 40 degrees, for the redshift intervals of 0.57 < z < 1.07 
(1500 < rc < 2500/i"^Mpc) and 1.07 < z < 1.78 (2500 < 
Tc < 3500 /i~^Mpc). For comparison, at rc = 3 /i~^Gpc, the 



simulation box-size subtends an angle of about 18 degrees 
on the sky. We find, as expected, that the intensity of the 
ISWRS and LAV maps drops with increasing redshift. By 
contrast, the intensity of the RS maps remains similar at 
the two epochs so that this becomes dominant effect at high 
redshift. 

We now consider the one and two point statistics mea- 
sured from the maps. The results are summarised in Figs. 1131 
- \17\ which correspond to projections over the following re- 
spective redshift ranges: < z < 0.17, 0.17 < z < 0.57, 
0.57 < z < 1.07, 1.07 < z < 1.78, and < ^ < 10.3. Start- 
ing with the bottom left-hand panels, we see that the his- 
tograms of the temperature fluctuations of the ISWRS and 
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Figure 13. Statistics of the ISWRS, LAV and RS maps shown at the top panel {0 < z < 0.17) of Fig. El Fig. [9] and Fig. [Tol Top-left: 
angular power spectra. The orange line at the top is the WMAP5 measurement of the CMB temperature power spectrum. The ISWRS 
power spectrum is shown in black line. The red line is the power spectrum of the map constructed by applying the linear approximation 
to the velocity field (LAV), i.e. equation The blue line is the linear theory prediction for the same redshift interval. Top-right: the 
ratio of the LAV and ISWRS in the top-left panel. Bottom-left: the histogram of pixel temperatures of the ISWRS (black) and LAV 
(red) maps. The number of pixels within each bin of 0.1 fiK has been divided by the total number of pixels. Bottom-right: the histogram 
of pixel temperatures of the residual map (RS). Ss is the skewness of the map temperature. 



LAV maps are highly non- Gaussian at low redshift (partic- 
ularly in Figs. [13] and [14]), but are closer to Gaussian at 
higher redshifts. At low redshift the distributions are sig- 
nificantly skewed and irregular due to the fact that only 
a small number of features contribute to these maps. It is 
to be expected that potential variations in our local uni- 
verse will also generate similar large-scale temperature fluc- 



tuations in the observed CMB and these too will be non- 
Gaussian. To quantify this effect on the observed CMB, it 
is necessary to investigate the contribution of the local envi- 
ron ment on scales up to a t least a few hundred Megaparsecs. 
(See iMaturi et alj (|2007l l for discussion of the contribution 
from within 110 Mpc of Earth.) At higher redshifts, the rel- 
atively small size of the simulation box means that the same 
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Figure 14. As Fig. 1131 but showing the statistics for the maps in the bottom panel of Figs [8] [9] and 1101 projecting over the range 
0.17 <z < 0.57. 



features will appear more than once in the maps, which will 
lead to some artificial residual non-Gaussianity, but only on 
large angular scales. 

Although the one point distributions are very similar 
there is a noticeable systematic difference between the his- 
tograms of the ISWRS and LAV effects at all redshifts. The 
ISWRS line is higher/lower than the LAV line in the nega- 
tive/positive tails of the distribution. Because the monopoles 
of the maps have been set to zero, so the means of the distri- 
butions are the same, these features are indicative of differ- 
ence in the skewness of the distributions. The bottom-right 
panel shows the one point distribution of RS effect, given 



by the difference map formed by subtracting the LAV from 
the ISWRS map. The resulting distribution, which has zero 
mean, is clearly skewed. Measuring the skewness of the RS 
distribution, 5*3 = Ms/c^^, where /is is the third moment 
about the mean and a is the standard deviation of pixel 
temperatures, we find 5*3 ^ 0.1 over all the redshift inter- 
vals. The positive skewness indicates a shift of the mode of 
the distribution to negative temperature values, consistent 
with our findings that the RS effect produces negative per- 
turbations in both overdense and underdense regions. This 
skewness is a further indication, already visible in the RS 
sky maps, that the RS effect is significantly non- Gaussian. 
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Figure 15. As Fig. 1131 but showing the statistics for the maps in Fig. [TT] projecting over the range 0.57 < z < 1.07, corresponding to 
the comoving distance from 1500 to 2500 Mpc. 



We show the results of the two point statistics, in the 
form of the angular power spectrum, in the top-left panels 
for all five redshift ranges in Figsll3l-I171 The expected (lin- 
ear) ISW angular power spectrum is shown as a blue line. 
Since the simulation we use to compute the maps is peri- 
odic, and there is no power on scales larger than the box, 
the angular power in our maps drops below the linear ISW 
prediction at the largest angular scales. The angular scale at 
which this departure occurs depends on the redshift range 
of the projection, and occurs at higher / with increasing red- 
shift. Our simulation box-size is 1000/i~^ Mpc on a side so 
there are waves in the simulation which are well described by 



linear theory. So, as expected, over intermediate values of /, 
we find excellent agreement between our ISWRS (and LAV) 
maps and the ISW prediction. We plot the ratio of the LAV 
and ISWRS power in the top-right panel. Using the LAV ap- 
proximation significantly underestimates the power at large 
scales. Comparing the LAV and ISWRS power spectra, we 
see, in Fig. [TS] that the LAV approximation underestimates 
the power by 20% at around / = 35. At higher redshifts the 
value of / at which the power is underestimated by 20% in- 
creases gently to around 90 fFig jlOp . For the full projection 
fFig jl7p a 20% deviation occurs at around / = 80. 

In Fig. [T7] we have also plotted the power spectrum 
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Figure 16. As Fig. [13] but showing the statistics for the maps in Fig. 1121 projecting over the range 1.07 < z < 1.78, corresponding to 
the comoving distance from 2500 to 3500 Mpc. 



of the kinetic Sunyaev-Zel'dovich (kSZ) effect taken from 
ICooravl (|200lh for a similar cosmological model. The ther- 
mal and kinetic SZ effects both arise from Compton scatter- 
ing of CMB photons off electrons in the ionized intracluster 
medium. The power spectrum of the thermal effect is about 
one order of magnitude great er than the kSZ effect at arc 
minute and degree scales fe.g. lCooravll200ll : lHu DodelsonI 
I2OO2I ). However, the transfer of the thermal energy of the 
electrons to the CMB photons produces a characteristic dis- 
tortion of the CMB spectrum and so, in multi-band CMB 
observations, the thermal SZ signal can be separated from 
that of the ISWRS. In contrast, the kSZ effect, which is due 



to the bulk motion of the ionized material, does not share 
this spectral signature and so cannot be separated from the 
ISWRS signal. At / ^ 80 where the RS effect begins to domi- 
nate over the linear ISW contribution, the kSZ effect is more 
than an order of magnitude weaker than the ISWRS effect. 
However, the kSZ effect has a much steeper power spectrum 
and becomes comparable to the ISWRS power at / ^ 225, 
corresponding to an angular scale of ^0.5 degree. Thus, on 
sub-degree scales, both the ISWRS and kSZ effects must 
be considered. We note that their spatial characteristics are 
quite distinct. A moving lump composed of both dark mat- 
ter and ionized gas moving along the line-of-sight towards 
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Figure 17. As Fig. 1131 but showing statistics of maps from z = 0.0 to z = 10.3, correspondin g to the comoving distance from to 
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Mpc. The cyan hue, in the top left panel, indicates the predicted kinetic SZ effect from lCooravl (|2001.) 



the observer will produce a kSZ hot spot while one moving 
away from the observer will produce a coldspot, but neither 
case will produce an RS distortion. In contrast, if the lump 
is moving in the transverse direction it will produce an RS 
dipole, but no kSZ distortion. 



5.1 Cold spots in the CMB? 

In this section, we compare the properties of cold spots 
produced by the ISWRS with the notable C MB cold spot 
found in the southern Galactic hemisphere (|Vielva et al.l 



Cruz et al.1 l2(305l; ^icEwen et al.1 l2005l: ICmz et al 



l2006|; iMcEwen et"alll2006l : ICruz et al.ll2007l : iMcEwen etH 
2008). In Sections 14.21 and we have illustrated how both 
convergent and divergent flows can generate negative tem- 
perature perturbations via the non- linear RS effect. For con- 
vergent flows into overdense regions this temperature decre- 
ment will be swamped by the increment caused by the linear 
ISW effect, except at high redshift where the linear effect 
diminishes more rapidly than the non- linear effect. The am- 
plitude of both the high- 2; decrements and low- 2; increments 
is expected to be of order 1 /xK or less (see, also the studies 
of non-linearity in cluster regions [Martinez- Gonzalez et all 
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Figure 18. Legendre transforms of the Spherical Mexican Hat 
Wavelet given by equation [8] for different widths of the filter scale, 
R. The filters are normalized so that 27r J ^^ijj-, R)d cos = 1 and 
^ is related to Wr(1) by equation (9). 
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gent flows both the linear ISW effect and the non-linear RS 
effect reinforce each other to produce decrements that have 
amplitudes of around 10 /iK at low redshift (e.g. Fig. [8]) with 
app roximately 10% co ntributed by the non-linear RS effect. 
ICruz et aP (|2005l ) estimate that the southern cold spot 



has a temperature decrement of 73 /iK and an angular ex- 
tent of 10° when smoothed by a Gaussian filter of width 
4 degrees. To suppress the contributions to this decrement 
from primordial perturbations of both smaller and larger an- 
gular scales and hence b e tter c haracteriz e the properties of 
the cold spot ICruz et al.1 (|2005l ) (see also lVielva et al.ll2004l : 
ICruz et al.1 [20061 . l2007l ) filter the sky map with a spherical 
Mexican hat wavelet (SMHW). The SMHW filter can be 
expressed as 



/2^N{R) 

where N(R) = (1 +_R^/2+_R^°-^ 



(8) 

2 tan(6 '/2) and is 



the polar angle ([Martinez- Gonzalez et al ] [2002l ). The effect 
of smoothing with this form of filter can be understood most 
easily in spherical harmonic space. In Fig jlSI we show the 
Legendre transform of the SMHW, 



7* 



Wr(1) = 2-k \ 'iiy,R)Pi{cos9)dcose, 



(9) 



wher e Pi is the Legendre polynomial of order / ([Page et al.[ 
[2003[ ) , for various smoothi ng scales, R. Note t h at these fil- 
ters are normalized as in [Vielva et al] ([2004[ ): [Cruz et al.[ 
([2005[ . [2006[ . [2007[ ) such that 27r / ^'^{y,R)dcosO = 1. We 
see that this set of filters retain sensitivity to modes in the 



range 5 < / < 40, while effectively suppressing the primor- 
dial C MB perturbations that peak at / = 200. [Cruz et aU 
(|2005[ ) find that when filtered in this way with an SMHW 
of = 4.17 degrees the temperature decrement of the cold 
spot in the WMAPl d ata is —16 mK[^, w ith similar results for 
= 4 and 5 degr ees (|Cruz et al.|[2006[ ) and with WMAP3 
([Cruz et al.[ [2007[ ). The probability of finding such a high 
deviation in Gau ssian si mulations of the CMB i s 1.85% 
([Cruz et al.|[2007[ ). but see [Zhang Hutererl (|2009[ ). There 
i s no evidence that the decrement is frequency dependent 
([Cruz et al.[ [2006[ , [2007[ ) and hence any contribution from 
the thermal SZ effect must be small. 

To compare our simulated sky maps with the proper- 
ties of this observed cold spot we first smoothed the ISWRS 
maps of Figs. 8, 11, and 12 using a Gaussian filter with the 
FWHM of 4 degrees. The resulting maps contain cold spots 
of ^ 10 degrees extent, but only with the amplitudes of 
up to —20 /xK, approximately a factor of three lower than 
the reported amplitude, —73 /iK, though, of course, the ob- 
served cold spot may contain contributions from primordial 
CMB perturbations, which are missing from out maps. We 
find that the majority of the deepest cold spots come from 
the redshift range 0.17 < z < 0.57, which corresponds to 
a range of comoving distance from the observer of 500 - 
1500 Mpc. The optimum smoothing scale is R — 6- 
10 degrees. 

A more stringent comparison with the observed cold 
spot can be made by applying the SMHW filtering to our 
ISWRS maps. An example with an R = 5 degrees SMHW 
is shown in Fig. [19] which shows the map resulting from the 
redshift range 0.17 < z < 0.57. It should be noted that the 
amplitude of the fluctuations in this map are much less than 
those in the G aussian filtered m aps reported above simply 
because of the ([Cruz et al.[[2006[ ) choice of normalization of 
the SMHW filter. We, again, find the typical size of coher- 
ent structures is comparable with that found in the obser- 
vations. However, the coldest region we find from the map 
after SMHW filtering is - -2 /xK, accounting for 12.5% of 
the reported CMB cold spot. The left-hand panels of Fig. [19] 
show the difference between the full ISWRS prediction and 
that of only ISW (LAV). Careful comparison shows that this 
non-linear RS contribution acts to enhance the depth of cold 
spots (an example of which is shown in the lower panels of 
Fig. I19p and to diminish the height of the hot spots. How- 
ever it is a modest effect, contributing only 0.1 /iK (5%) to 
the depth of the illustrated cold spot. 

Previous studies have mainly focused on the ISW(+RS) 
effect arising from a single void, i.e, a region that 
is significantly underdense relative to the mean. They 
have concluded that a void of size 100 Mpc at 
z < 1 would be needed to produce a cold spot 



tational effects ([Rudnick et al. 


2007: Inoue & Silk[ 2006[. 


2007[: [Tomita & Inoud [2008[: 


Masina & Notari" "2009a[|b[; 


Granett et al.[[2009[: [Francis & Peacock 2009: Granett et al. 


2009[) and that such a huge void 


is very unlikely to occur in 



the standard ACDM cosmology. Here we have generalised 



^ Often spectral filters are instead normalized such that their 
peak value Wji(lmax.) = 1. If this were done here, the temperature 



decrement would be ~ 120 /iK. 
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Figure 19. Full sky maps of the predicted CMB AT due to the ISWRS (top-left) and solely th e RS, i.e. ISWRS-L AV, (top-right) effect 
made from our simulation smoothed by a SMHW with R = 5.00 degrees, the same as Rg in ICruz et al The grid spacing is 

20 degrees in both longitude and latitude. The bottom panels show a zoomed in version of a 22 x 22 degrees patch taken from the top 
plots. These maps are made by ray-tracing through the simulation over the redshift interval 0.17 < z < 0.57, which corresponds to a 
range of comoving distance from the observer of 500 - 1500 Mpc. 



this comparison by simulating the correct 3D ACDM dis- 
tribution of evolving density perturbations and mimicked 
the observations by producing projected maps. However the 
conclusion remains the same that the ISW+RS effect in the 
standard ACDM cosmology does not produce cold spots as 
deep as the one discussed. 



6 CONCLUSION AND DISCUSSION 

Our aim in this paper is to investigate fully the gravita- 
tional effect that CMB photons suffer when passing through 
the evolving non- linear gravitational potential, of inter- 
vening large-scale structure. We have developed a method of 
using large cosmological N-body simulations to compute the 
time derivative of the potential, i>, along the past light-cone 



of an observer. By integrating along light rays we have cre- 
ated full sky maps, both using the full non-linear calculation 
of $ and an alternative in which the dark matter velocity 
field is assumed to be related to its density field by the nor- 
mal linear theory relation. By comparing the results of the 
two calculations we were able to assess fully the linear Inte- 
grated Sachs- Wolfe (ISW) and non-linear Rees-Sciama (RS) 
contributions to the induced CMB temperature fluctuations. 

In general, in a ACDM universe, the linear ISW effect is 
dominant at low redshift where the accelerating expansion, 
driven by dark energy, causes the decay of perturbations in 
the gravitational potential. The propagation of CMB pho- 
tons through these decaying potential wells and hills pro- 
duces hot and cold regions respectively with AT of the or- 
der of lO/xK on scales of hundreds of Megaparsecs. At low 
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redshift the non-linear (RS) effect produces only small scale 
perturbations to this large scale pattern. However, with in- 
creasing redshift the RS contribution decreases in amplitude 
much more slowly than does the contribution from the ISW 
effect, which vanishes as Q.m{z) approaches unity. Hence, 
the importance and scale of the RS effect beco mes larger 
at hig her r edshift, confirming the conclusions of ICai et al.l 
(I2OO9I ) and ISmith et all (l2009l ). We have shown that the 
origin of the RS effect is primarily the non-linear relation 
between the velocity and density field rather than the non- 
linearity of the density field itselfQ Our investigation of the 
RS effect has revealed three distinct non-linear phenomena 
that give rise to corresponding characteristic features in the 
temperature perturbation maps. 

• Dipoles are produced by the transverse motion of large 
lumps of dark matter, with typical sizes of tens of Mega- 
parsecs, much larger than the scale of individual galaxy clus- 
ters. 

• Convergent flows, on the scale of up to 100/i~^Mpc, 
around non-linear overdense regions give rise to cold spots 
of order /xK surrounded by hot filamentary shells. At high 
redshift these can be strong enough to dominate over the 
linear ISW effect and change the sign of the temperature 
perturbation centred on overdense regions. 

• Divergent flows around void regions are characterised 
by RS contributions consisting of hot rings around cold cen- 
tral regions, where the density contrast of the void is nearly 
saturated [5 ~ — !)• This is a small effect at low redshift, 
but acts to strengthen the cold spots produced by the linear 
ISW effect. 

Unfortunately, none of these effects can be easily detected in- 
dividually. At low redshift these phenomena make only 10% 
changes to the temperature perturbations predicted by the 
linear ISW effect. At very high redshift they are completely 
dominant, but their amplitudes are very low. We find that 
they contribute to cold spots of comparable physical scale 
to those reported in the literature, but their amplitudes are 
many times smaller. 

It may be possible to detect these large scale fea- 
tures induced by the non-linear velocity field by employing 
stacking techniques. The detect ability of RS kinematic fea- 
tures produced by nierging clus ters has been discussed by 
iRubifio-Martm et al.1 (|2004l ) and lMaturi et al.1 (|2007l ). They 
conclude that around a thousand clusters would be needed 
to detect the RS signal above the contaminating effects of 
the primordial CMB temperature fluctuations and instru- 
ment noise. 

However, we have found imprints of the RS perturba- 
tions on scales of a few tens of Megaparsecs, much larger 
than the merging cluster scale, and with slightly larger am- 
plitudes. Thus, these large-scale features might be more eas- 
ily detected, requiring the co-addition of fewer objects. 

From our all sky maps we find that the RS contributions 
to the overall power spectrum of temperature perturbations 



^ The fact that to model ^ and the growth of large-scale struc- 
ture requires more accurate modelling of the large-scale velocity 
field than is provided by linear theory may suggest that using 
redshift space distortions to measure the growth rate of density 
perturbations, ^, may also not be sufficiently accurate. 



become important for / > 80 (a few degrees) and completely 
dominate for / > 200. At still smaller scales the kinetic 
SZ effect is expected to dominate and at such scales a full 
treatment would have to incorporate additional modelling 
of this contribution. The RS contribution to the tempera- 
ture maps is strongly non- Gaussian with a skewed one-point 
distribution. In future work it will be interesting to inves- 
tigate the RS contribution to higher order statistics as its 
non-Gaussian characteristics might limit the ability to de- 
tect primordial non-Gaussianity in the underlying primary 
CMB fluctuations. Combining our full-sky maps with mock 
galaxy catalogues built from the same N-body simulations 
will be a powerful tool for developing cross-correlation tech- 
niques aimed at extracting the ISWRS signal from redshift 
surveys. 
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